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We have carried out a systematic study on the effect of Cu doping on nuclear, magnetic, and dielectric 
properties in Mni_vCutW 04 for 0 < x < 0.19 by a synergic use of different techniques, viz, heat capacity, 
magnetization, dielectric, and neutron powder diffraction measurements. Via heat capacity and magnetization 
measurements we show that with increasing Cu concentration magnetic frustration decreases, which leads to the 
stabilization of commensurate magnetic ordering. This was further verified by temperature-dependent unit cell 
volume changes derived from neutron diffraction measurements which was modeled by the Griineisen approx¬ 
imation. Dielectric measurements show a low temperature phase transition below about 9-10 K. Further more, 
magnetic refinements reveal no changes below this transition indicating a possible spin-flop transition which 
is unique to the Cu doped system. From these combined studies we have constructed a magnetoelectric phase 
diagram of this compound. 


PACS numbers: 75.85.+t, 61.05.F-, 71.27.-ta, 75.30.Kz 


I. INTRODUCTION 


A typical feature of multiferroic materials undergoing a 
transition to an elliptic spiral ferroelectric phase, is the exis¬ 
tence of spectacular magnetoelectric effects, such as the po¬ 
larization flops observed in TbMnOa IH and orthorhombic 
DyMn03 ||2l or the sign reversal of electric polarization which 
is revealed under magnetic field in TbMn205 l^]. The orienta¬ 
tion of the applied magnetic field with respect to the magnetic 
spins influences the stability range of the spiral phase and the 
electric polarization-flop process. This property was recently 
illustrated by remarkable magnetic field induced effects ob¬ 
served in ferroelectric phase of manganese tungstate MnW04 
in which applied held induces a polarization flop transition 10- 

S. 

In most of the recently discovered multiferroics, the 
ferroelectric polarization can be explained by the inverse 
Dzyaloshinski-Moriya effect ifsl- fl^ . where the induced elec¬ 
tric polarization of a single pair of spins 5’, , ^j separated by a 
distance vector fi j is given by llsll 


PpE ocrijx X . 


( 1 ) 


The required helical magnetic structure may arise from 
strong frustration. Since in addition the interaction, equa- 
tion[T] is only a second order effect, the ferroelectric polariza¬ 
tion is rather small in these materials. In the RMnOs iiiini] 
(R=rare earth) series and in MnW04 12, 3 the electric 
polarization is about two to three orders of magnitude smaller 


than in a classical ferroelectric perovskite such as BaTiOs. As 
a consequence the observation of electric-fleld-induced ef¬ 
fects in the magnetically ordered state is more difficult. Nev¬ 
ertheless, it was shown that it is possible in these chiral multi- 
ferroics to switch the magnetic order ^ the application of an 
electric field at constant temperature lll4l4l6ll . 


The crystal structure of MnW04 is monoclinic with space 
group P2 /c, made up of MnOe octahedra with high-spin Mn^^ 
(ds) ions and WOe octahedra with diamagnetic (do) 
ions Ini]. Recently it was found that MnW04 exhibits multi- 
ferroicity in which magnetism causes ferroelectricity, imply¬ 
ing a strong coupling between the two lll[il[il. MnW04 is 
one of the prototypical multiferroic materials exhibiting spin- 
current ferroelectricity J^]. It possesses a complex phase di¬ 
agram with 3 antiferromagnetic phases below 14 K namely 
AFl, AF2 and AF3 at zero magnetic held. AF2 is a ferro¬ 
electric (FE) phase, in which the net polarization is along the 
b axis which can be flipped to the a axis with the applica¬ 
tion of an external magnetic held. This is the first example of 
the ferroelectric polarization flop induced by magnetic fields 
in transition-metal oxide systems without rare-earth 4/ mo¬ 
ments. Taniguchi et al. showed that the stability of the mag¬ 
netoelectric domain walls in a canted magnetic held plays a 
key role in the directional control of the electric polarization 
flop phenomenon ||3l. From polarized neutron scattering mea¬ 
surements Sagayama et al. showed that an inverse effect of 
Dzyaloshinskii-Moriya interaction is the origin of the sponta¬ 
neous electric polarization in the spiral phase of MnW04 |@1. 
From superspace symmetry formalism it was shown that in the 
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AF3 phase, the modulations of two Mn atoms within the unit 
cell can have a cycloidal component with equal and opposite 
chiralities canceling their effects and hence no electric polar¬ 
ization is induced. Whereas in the AF2 phase, an additional 
second magnetic mode with the spin modulations breaks the 
symmetry relation between the two manganese atoms with 
chiralities of the same sign which add up to induce macro¬ 
scopic electric polarization il8 l. 

Recently it was found that the ferroelectric phase is com¬ 
pletely suppressed in MnW 04 by doping 10% iron on Mn 
site, which can be again restored with the application of 
a magnetic field. The absence of ferroelectricity (at zero 
field) in Mno, 9 Feo.iW 04 is explained by the increase of uni¬ 
axial spin anisotropy K 0. Evidence for the increase of K 
with Fe substitution was also derived from neutron scatter¬ 
ing experiments S. It was observed that in Mni^;t Mj:W 04 
(M=Mg, Zn and x < 0.3), the substitution of the nonmag¬ 
netic Mg“^ ions and Zn^'^ for the magnetic Mn^^ ions result 
in very similar effects on the magnetic and dielectric proper¬ 
ties of MnW 04 12in . These substitutions destabilized the non¬ 
polar magnetic structure AFl of MnW 04 but the AF3-to-AF2 
magnetoelectric phase transition was not affected. This indi¬ 
cated that the nonmagnetic dopant destroys neither the three- 
dimensional nature of magnetic interactions, nor the spin frus¬ 
tration within each || c-chain and between || c-chains along 
the a-direction. In this article we discuss the influence of 
doping of Cu ions on the nuclear and magnetic structure of 
MnW04. 


II. EXPERIMENTAL 


Polycrystalline powders of Mni_j^Cu;cW 04 (jc= 0.0-0.19) 
were prepared by conventional solid state route. Stoichio¬ 
metric amount of precursors, W20j,(99.9%), Mn02(99.9%) 
and CuO(99.99%) were ground well with a mortar and pes¬ 
tle, pressed into pellets and sintered in a furnace at 900 °C for 
12 hours in the presence of atmospheric air. This process is 
repeated to achieve homogeneous powder samples. All com¬ 
positions were confirmed to be phase pure from x-ray powder 
diffraction. Specific heat measurements were carried out on 
small pellets, using a physical property measurement system 
(Quantum Design) in the temperature range 3-300 K. Mag¬ 
netic measurements were carried out in a commercial physical 
property measurement system using vibrating sample magne¬ 
tometer option. To investigate nuclear and magnetic struc¬ 
ture, time-of-flight (TOF) neutron powder diffraction (NPD) 
was performed on 8 g of powder samples that were loaded 
in 8 mm diameter vanadium cans. Neutron data were col¬ 
lected at the Spallation Neutron Source (SNS) at Oak Ridge 
National Laboratory on the high resolution neutron powder 
diffractometer POWGEN |22]. Data were collected for the 
compositions of Cu, x=0.0, 0.05, 0.1 and 0.19 in the temper¬ 
ature range 1.5 - 300 K. Eor each temperature the data was 
collected using two different center wavelengths, 1.599 A and 
3.731 A or 4.797 A. The crystal and magnetic structure refine¬ 


ments were carried out from the NPD data using the Rietveld 
refinement program Full Prof i23 . Eor the dielectric mea¬ 
surements, opposite sides of pressed pellets were covered by 
silver paint, thus forming a parallel-plate capacitor. The mea¬ 
surements were done using an LCR meter (Agilent 4980). Eor 
cooling down to 5 K, a He bath cryostat was used. 


III. RESULTS AND DISCUSSION 


The thermodynamic signature of the transition between dif¬ 
ferent phases is usually detected by distinct anomalies in the 
temperature dependence of the heat capacity, (Cp). Multifer- 
roic materials with several subsequent transitions may show 
pronounced anomalies of Cp. In figure[T]we present the vari¬ 
ation of Cp with the temperature for Mni_^Cu;cW 04 . Eor ref¬ 
erence, the specific heat of a MnW 04 single crystal is also 
included. All the compositions show two anomalies at Tni 
and rN 2 - A third low temperature phase transition, Tns seen 
in the case of MnW(94 is already missing in the lowest Cu 
doped compound. This is associated with the phase transi¬ 
tion from the helical AE2 phase to the commensurate AEl 
phase 101. This result implies that with Cu doping a quick 
suppression of the AEl phase occurs, as a result ferroelectric 
AE2 phase is extended to the lowest temperature. Similar re¬ 
sults of quick suppression of the AEl phase were reported in 
Mni_;tCo,vW 04 and Mni_;tZn^W 04 l240 . 

Dielectric measurements were performed on compositions 
X = 0.05, 0.1, and 0.19. The temperature dependence of the 
dielectric constant (s'), normalized to the dielectric constant 
value at 5 K, is presented in figure|2l To exclude contributions 
from electrode polarization or grain boundaries, which can 
lead to so-called Maxwell-Wagner relaxations, here we show 
the results at a relatively high frequency of 105 kHz |2^ 271. 
Eor X = 0.05 and 0.1, several anomalies in e'(T) are found as 
indicated by the arrows in figures |2| a) and (b). Those around 



FIG. 1. (Colour online) The temperature dependence of specific heat 
Cp of Mni_i-CUvW 04 . Different curves are vertically offset by 8 units 
along Cp-axis, zero for each curve is defined by the horizontal dashed 
lines. Vertical dashed lines are guide to eyes indicating transitions 
Tni and 7 'n2. 
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FIG. 2. (Colour online) Temperature dependence of dielectric con¬ 
stant at 105 kHz, for Mni_,vCuvW 04 with: (a) x=0.05, (b) x=0.1 and 
(c) x=0.19. The arrows in (a) and (b) indicate the dielectric anoma¬ 
lies and arrows in (c) indicate the anomalies seen in specific heat 
measurement of the sample with x=0.19. 


12 K agree with the findings from specific heat (Fig.[T]i. With 
increased Cu concentration, these dielectric-constant anoma¬ 
lies become weaker. In addition to the two transitions ob¬ 
served from specific-heat measurements, a third transition is 
revealed by the dielectric measurements at Tx = 9 and 10 K 
for the x=0.05 and 0.1 compounds, respectively. From neu¬ 
tron diffraction measurements we later show that Tx is not as¬ 
sociated with the phase transition from AF2 to AFl phase as 
in parent MnW 04 . As indicated by the arrows in figure |2tc), 
the two broad shoulders in e'(T) found for x=0.19 seem to 
roughly agree with the specific-heat results but a clear deter¬ 
mination of transition temperatures from the dielectric exper¬ 
iments is not possible for this sample. 


TABLE I. Structural data for MnW 04 , Mno. 9 Cuo,iW 04 and 
Mno 8 iCuo.i 9 W 04 obtained from the NPD data collected at 300 K. 





Cu content (x) 




x = 0 

X = 0.1 

x = 0.19 

fl(A) 


4.8300(5) 

4.8120(4) 

4.7946(6) 

b(k) 


5.7597(6) 

5.7645(5) 

5.7694(7) 

c(A) 


4.9977(5) 

4.9838(4) 

4.9708(7) 

V(A^) 


139.009(3) 

138.208(2) 

137.453(3) 

m 


91.140(7) 

91.370(7) 

91.579(9) 

Atoms 

Mn/Cu 

ylb 

0.6861(4) 

0.6854(6) 

0.6875(1) 


B,so (A") 

0.5120(47) 

0.4970(57) 

0.429(11) 

W 

ylb 

0.1815(3) 

0.1806(36) 

0.1806(3) 


(A^) 

0.5070(38) 

0.5550(38) 

0.6620(47) 

Global weighted (%) 

5.73 

8.2 

8.4 


,0.045 - 


(D 

o 


E 0.040 




0.035 - 


I I I I I 


(a) 


-.-0.05' 
0.1 - 
-.-0.15- 
-t-0.19: 








>- 










10 15 20 25 30 35 

Temperature (K) 



Temperature (K) 



Cu content (x) 


FIG. 3. (Colour online) (a) Susceptibility calculated from the magne¬ 
tization data measured after a field cooling cycle with applied mag¬ 
netic field of 1 kOe. (b) Curie-Weiss fit (red lines) to the inverse sus¬ 
ceptibilities, a deviation from linear nature is seen below magnetic 
ordering temperature.(c) The frustration parameter as a function of 
composition. Inset shows the Curie-Weiss temperature as a function 
of composition obtained from the Curie-Weiss fits. Straight lines are 
guide to eyes. 


Magnetization measurements of all samples were per¬ 
formed under magnetic fields of 1 kOe. Thermal evolution 
of magnetic susceptibility of the samples at low temperature is 
presented in figurej^a). From the magnetic susceptibility data 
of samples x=0.05, 0.1 and 0.15 only one magnetic ordering 
temperature is discernible around 14 K and for sample x=0.19 
two anomalies are discernible. The thermal evolution of in¬ 
verse susceptibility obtained from the field-cooled magnetiza¬ 
tion were fitted with Curie-Weiss law as shown in figure[3b). 
Inverse susceptibility of Mni_xCuiW 04 follows Curie-Weiss 
law down to ~ 75 K, below which it deviates from the fitted 
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FIG. 4. (Colour online) (a) Temperature evolution of diffraction 
patterns of Mno. 95 Cuo.o 5 W 04 . Magnetic phase transition is dis¬ 
cernible with the appearance of additional incommensurate Bragg 
peaks around 13.5 K (horizontal dashed line), (b) Temperature evo¬ 
lution of part of the diffraction patterns of Mno.gCuo 1 WO 4 . Mag¬ 
netic phase transition is discernible with the appearance of addi¬ 
tional incommensurate Bragg peaks around 14 K. (c) Temperature 
dependence of diffraction patterns of Mno. 8 Cuo,i 9 W 04 . Two mag¬ 
netic phase transitions are discernible at ~17 K and ~11.5 K. 


curve and shows a marked deviation below ~ 15 K which cor¬ 
responds to Tni. The deviation of inverse susceptibility well 
above ordering temperatures indicates the presence of short- 
range spin fluctuations above Tn. Thermal evolution of Curie- 
Weiss temperature (0cw) and the frustration parameter calcu¬ 
lated as f - |0cwI/Tn is presented in figure|3c) as a function 
of composition. Indeed MnW 04 has been known to be a mod¬ 
erately spin frustrated system with the frustration parameter, 
/ = |0 cwI/Tn ~ 4.9, where 0cw i s appro ximately -71 K and 
the Neel temperature Tn = 13.5 K| 12,^. From figure|3c) it 
is clear that Cu doping increases the Curie-Weiss temperature 
and reduces frustration. 

The temperature evolution of diffraction patterns in a 
selected (7-space is presented in figure 0] A mag¬ 
netic phase transition is discernible based on the new in¬ 
commensurate Bragg peaks below 13.5 K and 14 K in 
Mn()g 5 Cuo.o 5 W 04 and Mno. 9 Cuo. 1 WO 4 , respectively. In the 
case of Mno. 8 iCuo.i 9 W 04 two transitions were observed 



FIG. 5. (Colour online) (a) Unit cell volume and monoclinic angle as 
a function of Cu concentration at temperature 300 K, lines are only 
guides, (b)-(d) The temperature dependence of unit cell volume for 
the Cu compositions x=0.05, 0.1 and 0.19. The curves are fits to 
equation |2] with the Griineisen approximation for the zero pressure 
equation. Insets in (b)-(d) show the deviation of unit cell volume 
from the fitted data at low temperatures. In the second inset of (d) 
the unit cell volumes are plotted together, to compare the magnitude 
of negative thermal expansion close to ordering temperature the y 
axes of x=0.1 and 0.19 are scaled by adding 0.375 and 1.155 
respectively to match the value of the sample x=0.05 at 100 K. The 
lines are guides to eyes. 


around 17 K and 11.5 K. The fundamental crystal structure 
of all compositions is monoclinic with space group Pljc. The 
unit cell volume and monoclinic angle at 300 K is presented 
in figure|5la), the refined structure parameters including errors 
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FIG. 6. (Colour online) Observed and calculated diffractions patterns and their difference for (a) Mno, 95 Cuo,o 5 W 04 at 1.5 K (b) Mno, 9 Cuo,iW 04 
at 1.5 K (c) Mno. 8 iCuo.i 9 W 04 at 1.5 K and (d) Mno. 8 iCuo,i 9 W 04 at 13.4 K. Circles are the measured intensities and the curve is the calculated 
pattern. Top and bottom vertical bars mark the positions of the nuclear and magnetic Bragg reflections respectively. Bottom curve is the 
difference between the measured and calculated patterns. 


are tabulated in table (U). With the increase in Cu concentra¬ 
tion, a decrease in the unit cell volume and increase in the 
monoclinic angle was observed. The volume change is in ac¬ 
cordance with the Vegard’s law, lattice volume decreases as 
X increases, because the ionic radius of Cu^^ is smaller than 
Mn2+ 0. The temperature dependence of unit cell volume 
is shown in figure |3b)-(d). All three compositions presented 
here show a negative thermal expansion, an increase in vol¬ 
ume with the decrease in temperature below 50 K. This effect 
seems to halt around 20 K and the volume starts to shrink be¬ 
low magnetic ordering temperature. The effect is more pro¬ 
nounced in lower Cu concentration and diminishes signifi¬ 
cantly with Cu doping as shown in the lower inset of fig¬ 
ure EJl. The anomalous behavior of the temperature variation 
of unit cell volume is due to the magneto-elastic effect associ¬ 
ated with the antiferromagnetic transition at the Neel temper¬ 
ature. In order to study the spontaneous magnetostriction it is 
necessary to determine the temperature variation of the lattice 
parameters and the unit cell volume in the absence of mag¬ 
netism. One way to determine the background temperature 
variation of the lattice parameter and unit cell volume is to 
extrapolate the paramagnetic high temperature region to low 
temperature by fitting with a polynomial function . This 
is only an approximation which works in some cases but in¬ 
volves some uncertainty. Alternatively, we used the Griineisen 
approximation for the zero pressure equation of state, in which 
the effects of thermal expansion are considered to be equiva¬ 


lent to elastic strain m. Within this approach the tempera¬ 
ture dependence of the volume can be described by, 

y(r) = yU(T)IBo + Vo (2) 

where, y is a Griineisen parameter, Bq is the bulk modulus and 
Vo is the volume at T = 0 K in the absence of magnetoelastic 
effect. By adopting the Griineisen approximation, the internal 
energy U{T) is given by 

UiT)-9Nk.T[-] I ^dv (3) 

where N is the number of atoms in the unit cell, is Boltz¬ 
mann’s constant and is the Debye temperature. By fit¬ 
ting the unit cell volume in the paramagnetic state we can 
get the physical parameters and Tq- From the present fit¬ 
ting procedure it was not possible to determine y and Bq but 
the ratio of y/Bo was set as variable. The result of the fit 
is shown as a green solid line in figure |5lb-d). Remarkably 
the fitted curves deviate from the experimental data at around 
50 K much above the magnetic ordering temperature, where 
the unit cell volume undergoes a negative thermal expansion. 
Below magnetic ordering temperature Tni the negative vol¬ 
ume effect seized and the unit cell volume is decreased with 
temperature. This is inferred as a clear indication of presence 
of magnetoelastic or magnetovolume effects in these system 
as well as underlying frustration, though the negative ther¬ 
mal expansion is significantly small for x=0.19 compound. 
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TABLEII. The little group of k=(a, 1/2, 7)=(-0.2183, 1/2, 0.476) 
in P2/c is Gk=Pc (with co-set representatives with respect to the 
translation group: l,c). The corresponding IRs are one-dimensional 
ri(l, c)=(l, e) and r 2 (l,c)=(l,-£), with -£=exp(7tiy). The basis 
vectors of the IRs of Gk are given below for the atoms Mn/Cu in the 
primitive unit cell numbered as: 1(1/2, v, 1/4) and 2(1/2, 1-y, 3/4) 
related by c-glide plane: c (x, -y-l-1, z-l-1/2), respectively. 


IR BV 

Atom 


m\\b 

BV components 

mile im\\a 


mi||e 

Fi 

1 

1 

0 

0 

0 

0 

0 


2 

-0.075 

0 

0 

0.997 

0 

0 


1 

0 

1 

0 

0 

0 

0 


2 

0 

0.075 

0 

0 

-0.997 

0 

<Pi 

1 

0 

0 

1 

0 

0 

0 


2 

0 

0 

-0.075 

0 

0 

0.997 

r2 fi 

1 

1 

0 

0 

0 

0 

0 


2 

0.075 

0 

0 

-0.997 

0 

0 


1 

0 

1 

0 

0 

0 

0 


2 

0 

-0.075 

0 

0 

0.997 

0 

fi 

1 

0 

0 

1 

0 

0 

0 


2 

0 

0 

0.075 

0 

0 

-0.997 


Temperature evolution of lattice parameters are similar to Co 
doped compound in which complex magnetic phase diagram 
is attributed to lattice changes 124 1. 


Representational analysis allows the determination of the 
symmetry-allowed magnetic structures that can result from 
a second-order magnetic phase transition, given the crys¬ 
tal structure before the transition and the magnetic propaga¬ 
tion vector (k) of the magnetic ordering. These calculations 
were carried out using the program Basireps included in the 
FullProf suite. First, the program iearc/r, also included in 
the FullProf suite, is used to determine the magnetic propa¬ 
gation vector at different temperatures. For x=0.05 and 0.1 the 
magnetic propagation vector was found to be k = {k^, kf) 
in whole temperature range. For x=0.19 the magnetic prop¬ 
agation vector in the temperature range 11.5-17 K was found 
to be k = (|,0,0) and below 11.5 K it is k = {k^, 5 , kf). 
While the magnetic propagation vector determines the modu¬ 
lation going from one unit cell to another, magnetic symme¬ 
try analysis is needed to determine the coupling between the 
symmetry related magnetic sites within one crystallographic 
unit cell. Basireps is used to compute all the allowed sym¬ 
metry couplings in the form of irreducible representations and 
their respective basis vectors. The Mn/Cu at the site 2f in 
the crystallographic space group P2/c, for the incommensu¬ 
rate magnetic propagation vector k = (k^, j, kf), is found to 
have two possible irreducible magnetic representations (IR) 
each having three basis vectors (BV). The computed results 
for the x=0.05 compound at 1.5 K with the propagation vector 
k = (-0.218, 0.476) are presented in table (I^. All possible 

combinations of the two allowed irreducible representations 
were tested against the measured data. Rietveld refinements 
clearly showed that only with the IR F 1 (with real and imagi¬ 
nary components) a successful! refinement of the data can be 
obtained. The propagation vector and the refined coefficient 


TABLE III. Refined unit cell parameters, magnetic propagation vec¬ 
tor and coefficients of basis vectors for Mni_j:Cuj:W 04 at 1.5 K. 


Cu content (x) 



X = 0.05 

X = 0.1 

x = 0.19 

Unit cell parameters 

fl(A) 

4.8145(6) 

4.8045(4) 

4.7856(8) 

b{k) 

5.7565(9) 

5.7583(6) 

5.7619(8) 

c(A) 

4.9860(8) 

4.9788(5) 

4.9653(9) 

pn 

91.20(1) 

91.30(9) 

91.49(6) 

Mn/Cu ylb 

0.6841(6) 

0.6840(6) 

0.6857(2) 

Biso (A^) 

0.274(6) 

0.184(2) 

0.203(1) 

Components of propagation vector (k,, = 

7 ) 


kx 

-0.218(3) 

-0.221(6) 

-0.223(4) 

K 

0.476(1) 

0.494(4) 

0.526(4) 

Coefficients of basis vectors 

Cl (real) 

3.41(1) 

3.28(1) 

2.79(9) 

C2(imaginary) 

-3.96(2) 

-4.08(9) 

-3.62(4) 

C3(real) 

2.95(1) 

2.91(1) 

2.32(3) 


of basis vectors. Cl, C2 and C3, for x=0.05, 0.1 and 0.19 at 
1.5 K is presented in table (HIH) . It should be noted that C1 and 
C3 are real coefficients while C2 is a pure imaginary coeffi¬ 
cient. For x=0.19 in the temperature range 11.5-17 K, with 
k = (^,0,0), four one-dimensional IRs were computed, for 
Fi and F 2 two BVs are allowed and for F 3 and F 4 only one 
BV is allowed. The Shubnikov groups (SG) of each IRs have 
the same symbol Pa2/c (in Belov-Neronova-Smirnova nota¬ 
tion), but they correspond to different magnetic Wyckoff posi¬ 
tions and origin shifts . The magnetic moments of the two 
atoms (l)(x, y, z) and ( 2 )(-x-i-l, -y-rl, -z-i- 1 ) in the paramag¬ 
netic unit cell have the following configurations for eanch rep¬ 
resentation: Fi:l(mr, 0 ,mj), 2 (-m^, 0 ,-m^); F 2 :l(mj^, 0 ,m^), 
F 3 : 1 ( 0 ,m,., 0 ), 2(0,-my,0) and F 4 :l( 0 ,m 3 „ 0 ), 
2(0,-m,., 0). Only the representation described by F 2 (SG- 
Pa2/c, Wyckoff site 4/ in the doubled unit cell) provides a 
successful! refinement of the data with 1.24(7) ps and 
OTj=1.28(8)//b for7’=13.5K. 

Refined NPD patterns are presented in figure | 6 ta-d). From 
the magnetic structure refinements Mno. 95 Cuo.o 5 W 04 was 
found to order at ~ 13.5 K, with the incommensurate propa¬ 
gation vector k = (-0.218, j, 0.477). The temperature depen¬ 
dence of components of incommensurate propagation vector 
k = (k^, i, kj) is presented in figure |7j a). With decrease in 
temperature a distinct change in k^ and ky is observed. In case 
of MnW 04 the magnetic phases AF2 and AF3 are incommen¬ 
surate with similar wave vector k = (-0.214, 5 , 0.457). Only 
AF2 phase with elliptical spin arrangement is ferroelectric 
which can be explained W spatial inversion symmetry break¬ 
ing spiral spin structure ISD . It should be noted that in the case 
of Mno. 95 Cuo.o 5 W 04 in the whole temperature range the struc¬ 
ture is found to be similar to that of AF2 in the parent com¬ 
pound. The magnetic structure in case of Mno. 9 Cuo. 1 WO 4 is 
similar to that of Mn 0 . 95 Cu 0.05 WO 4 but with a modified propa¬ 
gation vector k = (-0.224, 7 ,0.5) close to 14 K. The evolution 
of propagation vector with temperature in Mno. 9 Cuo. 1 WO 4 is 
presented in figure |7jb). Striking change in propagation vec- 
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FIG. 7. (Colour online) Temperature variation of components of in¬ 
commensurate propagation vector along x (k^) and z (fc) axes for 
(a) Mno. 95 Cuo.o 5 W 04 , (b) Mno. 9 Cuo. 1 WO 4 and (c) Mno. 8 iCuo.i 9 W 04 . 
Component of propagation vector along y is ky = 0.5. 


tor close to ordering temperature Tni in this case indicates 
that with increased Cu content of x=0.1 the propagation vec¬ 
tor along z-direction is nearly commensurate. From our pow¬ 
der diffraction measurements for x=0.05 and 0.1 compounds 
we don’t see any significant change associated with the tran¬ 
sition from Tni (AF3) to 7 ’n 2 (AF2) as seen from specific heat 
measurements. Considering very narrow temperature range 
between these two transitions it might be difficult to clarify 
this with our bulk powder measurements. We expect that for 
x=0.05-0.15, the magnetic ordering in the temperature range 
Tm<T<Tm should be collinear incommensurate phase as in 
MnW 04 151] with magnetic propagation vector similar to AF2 
phase. Further studies on single crystals with polarized neu¬ 
tron diffraction with smaller temperature steps will be a best 
tool for the detailed investigation of this structure. 

With further increase in doping in case of Mno. 8 iFeo.i 9 W 04 
we observe a commensurate magnetic (CM) structure at 
7’ni = 17 K with k = 0) which is similar to AF4 phase 

in parent compound. Below 11.5 K it undergoes another mag¬ 
netic phase transition to AF2 phase with propagation vector 
k = (-0.225, 7,0.531) which is modulated with temperature 
as shown in figure |7};. The incommensurate cycloidal (AF2) 


and the commensurate collinear (AF4) magnetic structure of 
Mno. 8 iCuo.i 9 W 04 is presented in figure Oa) and (b), respec¬ 
tively. The incommensurate structures for lower doping sys¬ 
tems is quite similar to the one presented in figure[8ja). 

From our comprehensive study of Mni_vCuvW 04 we are 
able to construct a tentative magnetic phase diagram as shown 
in figure |9] The magnetic phase diagram of Cu doped com¬ 
pound found to be much simpler than that of Co doped com¬ 
pound ll^[^ but very similar to Zn doped compound 1251. 
With higher doping concentration a collinear magnetic struc¬ 
ture is stabilized at higher temperatures. This is attributed to 
weakening of spin-phonon coupling and hence lower frustra¬ 
tion leading to a simpler magnetic ordering. From the neutron 
diffraction measurement it is clear that the low temperature 
phase below Tx which is observed from dielectric measure¬ 
ments is incommensurate cycloidal phase. Magnetic structure 
refinements confirmed that the magnetic phase below Tx (re¬ 
gion marked with gray lines in figure |9]l is not associated with 
the transition from cycloidal structure with magnetic vectors 
k = (kjc, 0.5, k^) to collinear structure with magnetic vector 
k = (0.25,0.5, -0.5) as seen in MnW 04 d. This leads to the 
inference that below Tx the magnetic structure undergoes a 
temperature induced spin flip transition with similar magnetic 
propagation vectors which is indistinguishable from powder 
diffraction measurements. The suppression of low tempera¬ 
ture collinear phase can be attributed to extremely sensitive 
exchange coupling and anisotropy constants with respect to 
perturbations 34]. In the present case chemical dop¬ 
ing plays the role of perturbations. In a recent report based 
on magnetization, specific heat and ferroelectric polarization 
measurements, Liang et al showed that by the substitution of 
lower spin (1/2) Cu^^ for the higher spin (5/2) Mn^^ ion the 
multiferroic phases of MnW 04 are strongly affected 11381] . The 
Cu substitution will introduce a low spin with different ex¬ 
change coupling and anisotropy constants affecting the mag¬ 
netic and ferroelectric states. This leads to the stabilization 
of ferroeletric spin spiral phase at low temperatures with in¬ 
creasing Cu content. The microscopic exchange interactions 
can be obtained through inelastic neutron scattering (INS) ex¬ 
periments investigating the magnetic excitations. According 
to INS experiments on MnW 04 , to explain the magnetic ex¬ 
citation spectrum, up to 11 different exchange pathways were 
required to fit the data provin g th e three dimensional char¬ 
acter of magnetic fluctuations 11351] . This three dimensional 
nature explains the robustness of cycloidal spiral order in Cu- 
doped MnW 04 , since percolation threshold for site dilution is 
much lower than for two-dimensional systems Based on 
a semiphenomenological Landau theory, authors in OTIl clar¬ 
ified the effect of different dopants on the phase diagram of 
Mni^;iM^W 04 (M=Fe,Zn,Mg). The origin of complex phase 
diagrams in these compounds is attributed to competition be¬ 
tween different superexchange interactions with contrasting 
behavior of doping with different ions. We expect that the 
theoretical analysis presented in 13711 should be compatible 
for Mni_^CuiW 04 as well. The temperature induced spin- 
reorientation remains to be unique to the present compound 
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FIG. 8. (Colour online) (a) ICM structure of Mno. 8 iCuo,i 9 W 04 with propagation vector k = (-0.225, ^,0.531), the number of visible unit cells 
along a, b and c direction are, 5, 2 and 20, respectively, (b) High temperature CM structure of Mno. 8 iCuo.i 9 W 04 with propagation vector (0.5, 
0, 0), two unit cells along all three axes are shown. Gray box indicate one unit cell. 


which requires further scrutiny. 


IV. CONCLUSIONS 

From our detailed investigation of the Mni_j:Cuj:W 04 for 
0 < X < 0.19 we have shown that substitution of Cu for 
Mn results in a reduction of the frustration. Also a reduc¬ 
tion in negative thermal expansion with the increased Cu dop¬ 
ing was observed which hints to a reduction in spin-phonon 
coupling with the higher Cu content. Temperature and dop¬ 
ing dependence of lattice parameters establish a strong depen¬ 
dence of magnetic structure on lattice changes. Both Tni and 
7 ’n2 increased with higher Cu content. This is in contrast to 



FIG. 9. Tentative phase diagram of Mni_vCU;tW 04 with phase 
boundaries obtained from specific heat (asterisk), NPD (inverted tri¬ 
angle) and dielectric (circle) measurements. The symbols in black, 
red and green corresponds to the transitions Tni, Tf ^2 and Tx, respec¬ 
tively. Blue asterisk indicates the transition from AF2 to AFl phase 
in MnW 04 . (PM-paramagnet, CM-commensurate magnet, ICM- 
incommensurate magnet, FE-ferroelectric and PE-paraelectric) 


Mni^iZnxW 04 1251]. The presence of third transition Tx is 
unique to the present compound. We note again, our NPD 
data confirms that the origin of Tx is not ICM to CM observed 
in MnW 04 at 7 'n3 . A possible origin of this phase transition 
is the temperature induced spin reorientation. Further single 
crystal neutron diffraction and electric polarization measure¬ 
ments are desirable to shed light on the nature of ferroelectric 
and magnetic ordering below Tx and in the region between 
Tni and T^i- 
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